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Abstract 

The large size limit of matrix integrals with quartic potential may be used 
to count alternating links and tangles. The removal of redundancies amounts 
to renormalizations of the potential. This extends into two directions: higher 
genus and the counting of "virtual" links and tangles; and the counting of 
"coloured" alternating links and tangles. We discuss the asymptotic behavior 
of the number of tangles as the number of crossings goes to infinity. 



1.1 Introduction and basic definitions 

This chapter is devoted to some enumeration problems in knot theory. For a 
general review of the subject, see [7]. Here we are interested in the application 
of matrix integral techniques. We start with basic definitions of knot theory. 

1.1.1 Knots, links and tangles 

We first recall the definitions of the knotted objects under consideration. A 
knot is a closed loop embedded in 3-dimensional space. A link is made of 
several entangled knots. An n-tangle is a knotted pattern with 2n open ends. 



2 



CHAPTER 1. KNOT THEORY - ZINN- JUSTIN AND ZUBER 



We shall be interested in particular in 2-tangles, where it is conventional to 
attach the four outgoing strands to the four cardinal points SE, SW, NW, NE. 

(5) (ID O >oo< JOT 

This figure depicts a knot, two links and three 2-tangles. 

All these objects are regarded as equivalent under isotopy i.e. under defor- 
mations in which strands do not cross one another, and (for tangles) open ends 
are maintained fixed. Our problem is to count topologically inequivalent knots, 
links and tangles. 

It is usual to represent knots etc by their planar projection with minimal 
number of over/under-crossings. There is an important 

Theorem (Reidemeister) . Two projections represent the same knot, link or 
tangle iff they may be transformed into one another by a sequence of Reide- 

v 
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Also, in the classification or the counting of knots etc, one tries 
to avoid redundancies by keeping only prime links. A link is 
non prime if cutting tranversely two strands may yield two dis- 
connected non trivial parts. Here is a non prime link: 




1.1.2 Alternating links and tangles 

We shall now restrict ourselves to the subclass of alternating knots, links and 
tangles, in which one meets alternatingly over- and under-crossings, when one 
follows any strand. 

For low numbers of crossings, all knots, links or tangles may be 
drawn in an alternating pattern, but for n > 8 (resp. 6) crossings, 
there are knots (links) which cannot be drawn in an alternating 
form. Here is an example of a 8-crossing non-alternating knot: 

One may show that asymptotically, the alternating links and knots are sub- 
dominant. Still the tabulation and counting of this subclass is an important 
task, as a preliminary step in the general classification program. 

A major result conjectured by Tait (1898) and proved in [16, 17] is the 

Theorem (Menasco, Thistlethwaite) . Two alternating reduced knots or links 
represent the same object iff they are related by a sequence of "flypes", 
where a flype is a combination of 
Reidemeister moves respecting the 
alternating character of tangles: 

We shall thus restrict ourselves to the (manageable) 





1 . 2. MATRIX INTEGRALS AND ALTERNATING LINKS AND TANGLES3 



Problem. Count alternating prime links and tangles. 

This problem was given a first substantial answer by Sundberg and Thistleth- 
waite in [22]. We will discuss in the rest of this text how the matrix integral 
approach has allowed to make significant progress building on their work. 

1.2 Matrix integrals and alternating links and tangles 
1.2.1 The basic integral 

Consider the integral over complex (non Hermitean) N x N matrices 

Zc = l dMe^i-' t + i (1.2.1) 

with dM = Yl^ j d^eMij dQmMij . It was proposed in the context of knot 
enumeration in [30]. 

According to the discussion of Chapter [?], its diagram- 
matic expansion involves oriented double-line propagators 
M 'j-: y 'kM*^ while its vertices may be drawn in a one-to- 
one correspondence with the previous link crossings, with, 
say, over-crossing associated with outgoing arrows. 

As usual, the perturbative (small g) expansion of the integral (1.2.1) or of 
the associated correlation functions involves only planar diagrams in the large 
N limit. Moreover the conservation of arrows implies that the diagrams are 
alternating: 

4\ ^ /\ /\ 

It follows from the discussion of Chapter XX that in the large limit 

limAr_).oo j^logZc = X] planar connected alternating i^^^^i , where |AutZ?| IS thc Order 

diagrams D with n vertices ' ' 

of the automorphism group of D. But going from complex matrices to hermitian 
matrices doesn't affect that "planar limit", up to a global factor 2. We thus 
conclude that, provided we remove redundancies including flypes, the counting 
of Feynman diagrams of the following integral over N x N Hermitean matrices 
M, for A^ oo, 




Z = |dMe^[-iFM^ + !l^1 (1.2.2) 



with dM = YlidMiiYl^^j d^eMij d^mMij , yields the counting of alternating 
links and tangles. 

1.2.2 Computing the integral 

The large A^ limit of the integral (1.2.2) may be computed by the saddle point 
method, by means of orthogonal polynomials or of the loop equations, as re- 
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viewed elsewhere in this book. 

In that N ^ oo hmit, the eigenvalues A form a continuous distribution with 
density u{X) of support [—2a, 2a], forming a deformed semi-circle law [3] 

= ^(1 - 2|a2 - |A2)V4a2-A2 (1.2.3) 



with a^ related to g and t by 



3^0^ - a^ + 1 = (1.2.4) 



and one finds that the large limit of the "free energy" F is 



p=i 



We recall that this formal power series of F, the "perturbative expansion of 
F" , is a generating function for the number of connected planar diagrams, (as 
usual, weighted by their inverse symmetry factor) 

OO ^ o<X) 




For future reference, we note that the asymptotic behavior of Fp as |j — )• oo is 

Fp ~ const (12)fp"'^/2 . (1.2.5) 

Also all the 2p-point functions ^(f M^p) = / p{X)X^P may be computed. 
We only give here two expressions that we need below, the 2-point function 

A = ^a2(4-a2)= —Q)— (1.2.6) 

and the connected 4-point function F 

1 \-/ 

T{g,t) = -^aHl-a^){2a'-5)= I |, (1.2.7) 

whose diagrams, after removal of redundancies, will count 2-tangles. The p-th 
term 7p in the g expansion of F behaves as 12Pp~^^'^. 

1.2.3 Removal of redundancies 

The removal of redundancies for the counting links and tangles will be done in 
two steps. First "nugatory" that are in fact irrelevant diagrams representing 
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patterns that may be unknotted, 



I , and "non-prime" diagrams 



^, may be both removed by adjusting t = t{g) in such a way that A = 
= 1. In the language of quantum field theory, this is a "wave function 
renormalisation" . 

We then find F{g) = F{g,t{g)) 



:i.2. 



24^. 



and V{g) = T{g,t{g)) = ^^^^g^ - 

In that way, one gets the correct counting of links up to 6 crossings and 
of 2-tangles up to 3 crossings. This is apparent on the following table where 
we have listed in (a) the first links with their traditional nomenclature; (b) 
the corresponding Feynman diagrams with their symmetry weight; (c) the first 
2-tangles in Feynman diagram notation. It appears that the diagrams of each 
of the last two pairs in (c) are flype equivalent. 

dg) 




i 



3, 



4i 



5^ 



(b) 



(D ^ ^ 





(') X >o< ^ >oo< X y A Ox >0 



For links the first flype equivalence 
occurs at order 6: 





The asymptotic behaviour Fp ~ const (27/4)^ exhibits the same "crit- 
ical exponent" —7/2 as in (1.2.5) but an increased radius of convergence, as 
expected. Likewise 7p ~ const (27/4)^ 

In a second step we must take the quotient by the flype equivalence. Sund- 
berg and Thistlethwaite [22] proved that the flype equivalence can be dealt with 
by a suitable combinatorial analysis. The net result of their rigourous analysis 
is that the connected 4-point function F can be deduced from T{g) by a suitable 
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change of variable: this final computation has been rephrased in [2(3] where it 
is shown that it can be elegantly presented as a coupling constant renormalisa- 
tion g ^ go- In other words, start from N ^ (J^tM"^ 
before. Then compute T{gQ) and determine gQ{g) as the solution of 

2 



fM4), fixt = t(go) as 



90 = a 



-1 



:i.2.9) 



(i-<7)(i + r(^/o)). 

then the desired generating function is T[g) = T{gQ). 

To show this, we introduce H{g), the generating function of "horizontally- 
two-particle-irreducible" (H2PI) 2-tangle diagrams, i.e. of diagrams whose left 
part cannot be separated from the right by cutting two lines. Its Feynman 
diagram expansion reads 



H 





Then the 4-point function F is a geometric series of H 

summing up to T=H/ (1 — H). 
Now under the flype equivalence 



Thus, with F, resp. 

H denoting generating functions of flype equivalence classes of prime tangles, 
resp. of H2PI tangles and if H' is the non-trivial part of H, H = g + H', F 
satisfies a simple recursive equation 

H' 




r = g + gr + 



l-H> 



(1.2.10) 



both relations being depicted as 



\ / 




Consider now the perturbative expansion of F((7o) computed for a new value 
g^ of the coupling constant, depicted as an open circle 



r(5o 
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If we want to identify it to T{g), it is suggested to determine = go{g) by 
demanding that go = S — "^gH' — . . . 



so as to remove the first flype redundancies, and the remarkable point is that the 
elhpsis may be omitted and that no further term is required. Indeed ehminating 
H' between the two relations (1.2.10) and go = g — "^gH' gives 



go = g( - 1 + = (1.2.11) 

" (i-5)(r + i)^ ^ ' 

which is equivalent to (1.2.9) and also to relations found in [22]. In the case of 
the matrix integral (1.2.2), it is convenient to parametrize things in terms of 
A = . ^ i . One finds 

~ (A-2)U-A) 

r = ^ (1.2.12) 

4 

4(A-2) 

go = — ^3 , (1.2.13) 

where F is the wanted generating function of the number of flype-equivalence 
classes of prime alternating 2-tangles. Eliminating T and go between the three 
latter equations results in a degree five equation for A 

AA^ io^ —20 — 1) 

A^g - 6A'^g + ^ - 2,2A^ + 64A - 32 = (1.2.14) 

5-1 

of which we have to find the solution which goes to 2 as 5 — )• 

A = 2 + 2g + Qg'^ + 20g^ + 78/ + 334/ + 1532/ + 73725'^ + 36734/ + 187902/ + 

This then gives for F the following expansion (given up to order 50 in [22]) 

f{g) = g+2/+4/ + 10/+29/ + 98/+372/ + 1538/ + 6755/ + --- (1.2.15) 

and the asymptotic behaviour of the p-th. order of that expansion reads 

/ 101 + V2100l\^ ./2 
7p~ const ( ^ j (1.2.16) 

with again the same exponent —5/2 but a still increased radius of convergence. 

At this stage, we have merely reproduced the results of [22]. Our matrix 
integral approach has however two merits. It simplifies the combinatorics and 
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recasts the quotient by flype equivalence in the (physically) appealing language 
of renormalization. For example using the results of [3], one may easily compute 
the connected 2£-function which counts the number of flype-equivalence classes 
of prime alternating ^-tangles [26] 

T,, = ^{A-2Y-\3i-2-{i-l)A) 



ce+i - 



^ y (-4)^ 
'+1 ^ ^ ' 



and the numbers of 3- and 4-tangles up to 9 crossings are given by 

Tg = 3^2 + iV + 51^^ + 1865^ + 708^^ + 2850g^ + 12099c/*^ + 53756c/9 + • • • 

(1.2.17) 

Tg = 12g^ + 90^"^ + 468/ + 2196c/^ + 10044g'^ + 46170/ + 215832/ + • • • 

(1.2.18) 

Our approach also opens the route to generalizations in two directions: 

• higher genus surfaces and "virtual" links. 

• counting of "coloured" links, with a potential access to the still open 
problem of disentangling knots from links. 

This is what we explore in the next two sections. 



1.3 Virtual knots 
1.3.1 Definition 

The large N "planar" limit of the matrix integral (1.2.1) has been shown to 
be directly related to the counting of links and tangles. It is thus a natural 
question to wonder what the subleading terms in the N~~'^ expansion of that 
integral, i.e. its higher genus contributions, correspond to from the knot the- 
oretic standpoint. If one realizes that ordinary links and knots may always be 
deformed to live in a spherical shell 5^ x /, where the interval / is homeomor- 
phic to [0,1], one is ready to see that higher genus analogues exist. In fact, 
these objects may be defined in two alternative ways. 

First, as just suggested, they are curves embedded in a "thickened" Riemann 
surface S := S x [0, 1], modulo isotopy in S, and modulo orientation-preserving 
homeomorphisms of S, and modulo addition or subtraction of empty handles. 

But one may also focus on the planar representations of these objects. This 
leads to the concept of virtual knot diagrams [12, 13]. In addition to the 
ordinary under- and over-crossings, one must introduce a new type of virtual 
crossing, which somehow represents the crossing of two different strands that 
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belong to different sides of the surface but are seen as crossing in the planar 

projection. Thus virtual knots diagrams are made of ^K^, and 

virtual links and knots are equivalence classes of such diagrams with respect to 
the following generalized Reidemeister moves 




That the two definitions are equivalent was proved in [4, 15]. See [8] for a table 
of virtual knots. 

Virtual alternating links and tangles are defined in the same way as in 
section 2: along each strand, one encounters alternatingly over- and under- 
crossings, paying no attention to possible virtual crossings. 



Here is a virtual link depicted in several alternative ways: 




in (a), using ordinary and virtual crossings; in (b), three equivalent represen- 
tations on a Riemann surface. As illustrated by this example, in the thickened 
Riemann surface picture, the counting should be done irrespective of the choice 
of homology basis or of the embedding of the link/knot. But this is precisely 
what higher genus Feynman diagrams of the matrix integral do for us! 

Remark: there is a notion of genus for knots (minimal genus of a Seifert 
surface) which is unrelated to the genus defined above (genus of the surface S). 
For the former notion in the context of knot enumeration, see [20, 21]. 

1.3.2 Higher genus contributions to integral (1.2.1) 

We thus return to the integral (1.2.1) over complex matrices 

Z{g,t,N) = I dMe^H f MMt + | f (MMt)^] 
and compute F{g,t,N) = -^log Z{g,t, N)/Z{0,t, N) in an N""^ expansion 

oo 

F(5,t,iV) = ^iV2-2"FW(5,t) 

h=0 

F^^\g,t) receives contributions from Feynman diagrams of genus h. F^^^ is (up 
to a factor 2) what was called F in the previous section. 

i7(l) 

was computed 

in [18], F^'^^ and F^^^ in [2] and [1]. From F one derives the expressions of 
A = I — ^ and F = 2^ — 2A^. Moreover the first two terms in the g power 
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series expansion of any F^^\g) are easy to get [31] and provide some additional 
information. 

As before, we remove the non prime diagrams by imposing that A{g, t{g, N),N) 
= 1, which determines t = t{g, N) as a double g and 1/N'^ expansion. One then 
finds the generating function of prime 2-tangles of minimal genus h, r^^\g), as 
the Af2-2/i iQj.yj^ jj-^ ii^Q 1/Af2 expansion of T{g,N): 

r(0) (g) = g+2 g2+6 g^+gi g5_,_408 g6_,_i938 g7+gQU g8 +49335 g9+260130 

r(i) (g) = g+8 g2+59 g3+420 g*+2940 9^+20384 g'5+140479 /+964184 +6598481 +45059872 g'^°+--- 
r(2)(g) = 17g3+456 g4+7728 g5+i04762 gS+i240518 9^+13406796 g8+135637190 g9+1305368592 g'^° + --- 
r(3) (g) = 1259 9^+62072 gS+1740158 g^+36316872 3^+627368680 g»+9484251920 g'^° + --- 
rW(g) = 200589 g'^+14910216 9^+600547192 33 + 17347802824 c/io+--- 
r<-^) (g) = 54766516 3^+5554165536 g^°+--- 



1.3.3 Table of genus 1, 2 and 3 virtual links with 4 crossings 




Table of virtual knots and links with 4 crossings. Objects are not distin- 
guished from their mirror images, see [31] for details. 

1.3.4 Removing the flype redundancies. 

The first occurences of flype equivalences occur in tangles with 3 crossings: 

ix-XiC X(r><-x>ax 

It has been suggested [31] that it is (necessary and) sufficient to take the quo- 
tient by planar flypes, thus to perform the same renormalization g — )• go{g) as 
for genus 0. In other words, we have the 
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Generalized flype conjecture: For a given (minimal ) genus h, T^^^ (g) = 
T^'^^go) is the generating function of Eype-equivalence classes of virtual alter- 
nating tangles. 

Then denoting by T^'^^g) = T^'^^go) the generating function of the number 
of flype equivalence classes of prime virtual alternating 2-tangles of minimal 
genus h, r(0)(5) is what was called r{g) in section 2, Eq. (1.2.10), while 

r<^) (g) = g+S g2+57 g''+2512 g^+W158 g'*+102837 g^+649862 g» +4086137 g'J+25597900 g^° + --- 

r<^) (s) = 17 g=*+456 g'^+7626 g^+loOQlO 9^ + 1155636 g''+11987082 9^+115664638 g9+1056131412 gi"+--- 
r^^Us) = 1259 9^+62072 g^+1727568 g''+35546828 g**+601504150 9^+8854470134 g'^°+--- 
fW(^g) = 200589 /+14910216 g*+597738946 9^+17103622876 5i''+--- 
f{5)(g) = 54766516 3^+5554165536 c/io+--- 



The asymptotic behavior of the number of inequivalent tangles of order p is 



^ , 101 + V2l00T \ ^.(,_,) 



p 

5 / 



ip 



40 



In [31], this generalized flype conjecture was tested up to 4 crossings for 
links and 5 crossings for tangles by computing as many distinct invariants of 
virtual links as possible. We refer the reader to that reference for a detailed 
discussion. No counterexamples were found. 



1.4 Coloured links 

1.4.1 The bare matrix model 

Let us first describe the "bare" model that describes coloured link diagrams. 
Since we are only interested in the dominant order as the size of the matrices 
N goes to infinity, we can consider, as was argued in section 1.2.1, a model of 
Hermitean matrices (as opposed to the complex matrices that were necessary 
in section 1.3 for virtual tangles). 

Let us fix a positive integer r - the number of colours - and define the 
following measure on the space of r Hermitean matrices Ma ■ 

n dM, exp iVtr ( - - + | ^ (MaM^f) (1.4.1) 

a=l \ a=l a,b=l j 

This measure has an 0(t) symmetry where the matrices are in the funda- 
mental representation of 0(r). 
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Figure 1.1: The two types of correlation functions of the 0(t) matrix model. 



Expansion in perturbation series of the constant g produces 
the following Feynman diagrams: they are fat graphs (planar 
maps) with vertices of valence 4, in which the colours cross 
each other at each vertex, see the figure. The summation over 
0(r) indices produces a factor of r for every colour loop. 
Thus, we have the following double expansion in g and r: 

lim ^'^^ ^ \ ^ ^ ^number of vertices(D) ^number of loops(D) 




AT^oo Af2 ^ \AntD 

4-valent diagrams D 



g ^ ' r 



where it is understood that the number of loops is computed by considering 
that colour loops cross each other at vertices. In other words, the model of 
coloured links gives us more information than the one-matrix model because 
it allows for a "refined" enumeration in which one distinguishes the number of 
components of the underlying link. 

Note that F is at each order in g' a polynomial in r, so that we can formally 
continue it to arbitrary non-integer values of r. 



Observables 

Let be the set of pairings of 2k points (sometimes called "link patterns"), 
that is involutions of {1, . . . , 2k} without fixed points. 

To each given link pattern vr of 2k points one can associate the quantity /jr 
(/ stands for "internal connectivity") as follows. It is the generating series of 
the number of alternating 2A;-tangle diagrams (or simply, of 4-valent fat graphs 
with 2k external legs) with a weight of r per closed loop and a weight of g 
per vertex, in such a way that the connectivity of the external legs, which are 
numbered say clockwise from 1 to 2k, is represented by vr (assuming as usual 
that colours cross at each vertex). See Fig. 1.1(a). 

From the point of view of the matrix model these observables 1,^ are not so 
natural. In principle one can define them as follows: 

^ ( 77 ^^"^ ■ ■ ■ ^"^'^ / ai = aj <^ j = i or j = 7r(i) 
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By 0(r) symmetry, the result is independent of the choice of the Oj as long as 
they satisfy the condition above, i.e. that indices occur exactly twice according 
to the link pattern vr. However this formula only makes sense if r > k. 

A more natural quantity in the matrix model is the "external connectivity" 
correlation function E.,^, which is defined in a very similar way: 

The only difference is that this time one sums over all Oj (which might produce 
additional coincidences of indices, and in fact always will if r < A;). 

The graphical meaning of E^^ is that it is the generating function of tangle 
diagrams with 2k external legs and prescribed connectivity outside the diagram, 
cf Fig. 1.1(b). Closing the external legs will produce closed loops which must 
be given a weight of r. However, crossings outside the diagram should not be 
given a weight of g. 

Noting that all the diagrams that contribute to E.,^ must have a certain 
internal connectivity, we can write 

E-K = ^ GtT.Tt'-^Tt' (1.4.2) 

tt' 

The coefficients G^^jr' are nothing but the natural scalar product on link pat- 
terns of same size 2k, defined as follows: 

G^inumbcr of cycles of tt o tt' „ „t ^ tj 

7r,7r' = ^ 7r,7r G ^2fc 

Graphically, it corresponds to gluing together the two pairings and giving a 
weight of r to each closed loop that has been produced. 

As a consequence of the formulae presented below, for positive integer r 
and k > T, the matrix G has zero determinant and formula (1.4.2) cannot be 
inverted in the sense that the E.,^ are actually linearly dependent. For example 
at T = 1 there is really only one observable per k (with one colour one cannot 
distinguish connectivities). It is however convenient to introduce the pseudo- 
inverse W of G, that is the matrix that satisfies WGW = W and GWG = G. 
The definition of G still makes sense for non-integer r, in which case G is 
invertible and W = G~^. We now sketch the computation of W following [5] 
(where it is called the Weingarten matrix, in reference to [2.')]). See also [-32] for 
a recursive way to compute W for generic r. 

G is a {2k — 1)!! x (2k — 1)!! symmetric matrix, with the property that it 
is invariant by the action of the symmetric group, where the latter acts on 
involutions by conjugation (cr • vr = avra"^): G^..,^^a--TT' = G-wy'^ or equivalently 
aG = Ga for all a G S2k- Furthermore, one easily finds that C[P2fc] contains 
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exactly once every irreducible representation of associated to a Young dia- 
gram with even lengths of rows. Thus, G is a linear combination of projectors 
onto these irreducible subrepresentations, which are of the form 

where A is a Young diagram with 2k boxes (A = 2fi for the projector to be 
non-zero) and is the associated character of the symmetric group. 

Finally one can write G = c^-P^^ where ^ is a Young diagram with k 
boxes, and the coefficients can be computed [5, 28]: 

Cm= H {T + 2j (1.4.3) 

Therefore, the pseudo-inverse W of G can be written as 

W= c~^P^^ (1.4.4) 



Loop equations 

Loop equations are simply recursion relations satisfied by the correlation func- 
tions of our matrix model. They can in fact be derived graphically without any 
reference to the matrix model, in which case the parameter r can be taken to be 
arbitrary (not necessarily a positive integer). We recall that we limit ourselves 
to the dominant order as — )■ oo. 

The recursion satisfied by E^^ is illustrated on Fig. 1.2. Start with one of the 
external legs (say leg numbered one), and look at what happens to it once one 
moves inside the "blob" . There are two possibilities: (i) it reaches a crossing, in 
which case one gets a factor of g and a new correlation function Et^' where vr' is 
obtained from vr by adding one arch around the leg number one; or (ii) it goes 
out directly and connects to the external leg 2i, i = 1, . . . , k (possibly creating 
a loop and therefore a factor of r if 7r(l) = 2i). This second situation is more 
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complex because naively the two blobs created by cutting the initial blob into 
two may still be connected by say 2ii lines. Let us consider the two limiting 
cases. If £j = we simply have two disconnected blobs and the contribution is 
Et^^Etj^ where tti and 7r2 are connectivities of size 2(z — 1) and 2{k — i). On the 
contrary if = k — 1, the two blobs are fully connected to each other according 
to a certain permutation a G 'S2(^k-i) it is clear that internal connectivity 
for one becomes external connectivity for the other, so that the contribution is 
of the form X^^^gp^ ^. ^ ^-Kilcr-jn- The crucial remark is that one can rewrite 
this as 7r2eP2(fe i) ^cr-7ri,7r2-^7ri-£'7r2 sveu if G is non-invertible. Indeed the 
Etj, due to formula (1.4.2), live in the image of G and therefore one can ignore 
the zero modes of G {G being symmetric, its image and kernel are orthogonal). 
In the general case in which there are 2£i connections between the two blobs 
with associated permutation a £ ^2^-, one has to break these connections by 
using the matrix W for link patterns of size 2£i ; calling vri (pi ) the connectivity 
of the first blob in which the 2ii legs connecting it to the other blob have been 
replaced with the link pattern ■ pi of size 2ii, and similarly for 7r2(p2) and 
link pattern p2 , we get an expression of the form 

k 

E^ = gE^, +Y,r'<''''' E ^pi,P2^.i(pi)^.2(p2) (1-4.5) 

(where VF0 = £'0 = 1). This equation allows to calculate the E.,^ iteratively, 
in the sense that to compute the l.h.s. at a given order, the E-,^ appearing in 
the r.h.s. are either needed at a lower order in or at the same order in g but 
have fewer external legs, than the E.,^ in the l.h.s. 

1.4.2 Removal of redundancies and renormalized model 

As in section 1.2.3, we now discuss how to go from the counting of (coloured) 
alternating link diagrams to the counting of actual (coloured) alternating links, 
that is up to topological equivalences. We recall that the process involves two 
steps: removal of nugatory crossings and consideration of prime tangles only, 
which amounts to a wave function renormalization (i.e. renormalization of the 
quadratic term of the action); and inclusion of flypes, which amounts to a 
renormalization of the quartic term of the action. However, a crucial difference 
with the model discussed in section 1.2.3 is that in the 0(r) model of coloured 
links one can introduce not just one, but two 0(r )-invariant quartic terms: 
besides the already present term of the form tr j(MaM{,)^, one can also 
have another term of the form tr ^^^M^M^, and one expects that this term 
will be generated by the renormalization [25] . We now summarize the equations 
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that we find. We start from the measure 



T 



n dM^ exp ( AT tr ( - ^ + I (M.M,)^ + | J] M^M^ 

a=l \ a=l a,b=l a, 6=1 

(1.4.6) 

The Feynman rules of this model now allow loops of different colours to "avoid" 
each other, which one can imagine as tangencies. The loop equations of this 
model generalize in an obvious those of section xxx and will not be written 
here. 

Next we define the following correlation functions A and ro,±: 



A = ^(i2)/r 
To = (i5^(i2)(34) - T(r + 1)A^)/ 



n 



r± = -^(12)(34) - lb /(13)(24) 

The are not directly defined in the matrix model, but the E.^ are - in fact 
the 4-point functions E.^ are obtained by differentiating the free energy with 
respect to gi and 52- But from the formulae of section 1.4.1 one can check that 
the Weingarten function for 4-point functions is invertible for r 7^ 1, — 2 (special 
cases which only require Fq, as discussed in detail in [26], and which we exclude 
from now on). Thus one can deduce the 4-point /^r from the Ej^. Also define 
the auxiliary objects (generating series of horizontally two-particle irreducible 
diagrams) 

Hn = l ^ 



(i-5)(i + ro) 



(iT<7)(i + r±) 

Then the equations to impose on the bare parameters gi, g2 and t as func- 
tions of the renormalized coupling constant g are 

A = 1 (1.4.7) 

gi = g{l -H+- H^) (1.4.8) 
92 = -g{Ho/T + (1/2 - 1/t)H+ - (1.4.9) 

Up to order 8, we find 

gi =g-(2g4 + (2+2T)g5+(l4+2T)gS+(26+16T+2T2)g7 + (134+56T+2T2)g8)+... 

92 = -(9'^+S''+335+{5+2r)g6 + (27+5r)g7 + (89+32T+T2)g8)+- 
t = l+2s+Tg2_2Tg^-63''-(8+10T)g^-{38+16T+3T2)g6 

-(104+86r+14T2)g7_(4io+338T+56r2+2T3)g8+... 
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Composing these series with the correlation functions allows to produce 
generating series for the number of coloured (prime) alternating tangles with 
arbitrary connectivity. For example, we find for 4- and 6-tangles (we only 
mention one pairing per class of rotationally equivalent pairings): 



-^("12) (34) 


= 9^+9^ + (3+r)g'l+{9+T)g^+{21+llr+r2 




■{101+32t+t2)/ 


-'"as) (24) 


= g+2g3+2g''+(6+3r)g5 + (30+2r)g« + (62- 


f40T 


+2r2)g7 + ... 


''^(''14) (25) (36) 


= 2g3+l8g5+i8g64.(i56_,_24T)g^+--- 






''^(''14) (26) (35) 


= 9^+75''+6g5 + (53+8T)g« + (154+6r)/+- 






■''(''12) (35) (46) 


= 2g3+2g4+(i6+2r)g5+{42+2T)g'5+{1714 


-44t- 


h2r2)g7 + ... 


■''(''14) (23) (56) 


= 4g4+8g^ + (42+7T)/+{156+14T)/+--- 






■''(''12) (34) (56) 


= 3g''+9g5 + (41+7T)/+(168+21r)g^+--- 







The superscript c means we are considering the connected generating series 
(corresponding to tangles which cannot be broken into several disentangled 
pieces) e.g. /('i2){34) = Ai2){34) - 1, -^('i3)(24) = Ai3){24), etc. 

In particular, note that 2/^'^i2){34) fi3)(24) at r = 1 reproduces Eq. (1.2.15), 

and similarly ^(14) (25) (35) +3^(^14) (26) (35) +^-^("12) (35) (46) +^-^("14) (23) (56) +^-^(''12) (34) (56) 
at r = 1 reproduces Eq. (1.2.17). 



1.4.3 Case of r = 2 or the counting of oriented tangles and links 

In the case of the (renormalized) matrix model with r = 2 matrices, a subset 
of correlation functions can be computed exactly, including the two- and four- 
point functions which are necessary for our enumeration problem. Instead of 
giving two colours to each loop, one can equivalently give them two orientations: 
not only does this give a nice interpretation of the enumeration problem as the 
counting of oriented tangles, but it is also the first step towards the exact 
solution of the problem. Indeed, as shown in [29], this reduces it to the solution 
of the six-vertex model on dynamical random lattices, which was studied in 
[24, 14]. 

The explicit generating series are given in terms of elliptic Theta functions 
and will not be given here; even their asymptotic (large order) behavior is 
somewhat non-trivial to extract, and we quote here the result of [29]: if jp is 
the p*^ term of one of the four-point correlation functions, 

7p ~ est g^^p (logp) 

where Qc is the closest singularity to the origin of these generating series; 
l/S'c ~ 6.28329764. Though the latter number is non-universal, the subleading 
corrections are; they correspond to a c = 1 conformal field theory of a free 
boson coupled to quantum gravity. 
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1.4.4 Case of r = or the counting of knots 

A case of particular interest is the limit r — )■ of the matrix model (1.4.1). This 
can be considered as a "replica limit" where one sends the number of replicas 
to zero. Alternatively, the r — t- matrix model can be written explicitly us- 
ing a supersymmetric combination of usual (commuting) and of Grassmannian 
(anticommuting) variables, see [20]. 

The observables are defined as follows: 

E^ = lim{-E^) 

T^O T 

that is they correspond to tangles which, once closed from the outside, form 
exactly one loop (i.e. form knots as opposed to links). 
The loop equations of the bare model become 

,P2^1Ti{pi)Ett2[p2) (1.4.10) 

2=1, pl,p2&P2l, 
2i7^-n-(l) 

where W is the pseudo-inverse of G = \\mr^o{T~^G). Note that according to 
(1.4.3,1.4.4), the factor cancels the trivial zeroes of G at r = 0. These 
zeroes are simple for diagrams with A3 < 1; the remaining diagrams (in size 
n > 6) have higher zeroes, making G non-invertible. 

Though Eqs. (1.4.10) cannot be solved analytically, it is worth mentioning 
that they are easily amenable to an iterative solution by computer; in fact 
the resulting algorithm is notably better than the transfer matrix approach of 
[11, 10], and one finds for example for the two point function A = E(i2) the 
following power series in g: (using a PC with 8 Gb of memory and 24h of CPU) 

1,2,8,42,260,1796,13396,105706,870772,7420836,65004584,582521748,5320936416, 
49402687392,465189744448,4434492302426,42731740126228,415736458808868, 
4079436831493480,40338413922226212,401652846850965808,4024556509468827432, 
40558226664529024000,410887438338905738908,4182776248940752113344, 
42770152711524569532616,439143340987014152920384,4526179842103708969039296... 

The objects being counted by this formula are also known as self-intersecting 
plane curves or long curves, see e.g. [0]. 

One can similarly take the limit r — t- in the renormalized model. However 
little is known beyond the general facts mentioned above for arbitrary r. 

1.4.5 Asymptotics 

The most interesting unsolved question about the 0(r) matrix model of coloured 
links and tangles concerns the large order behavior of the generating series in 
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the coupling constant g, i.e. the asymptotic number of coloured alternating tan- 
gles as the number of crossing is sent to infinity. If one considers, in the spirit 
of chapter that the model represents a statistical model on random lattices, 
then it is expected that the model is critical for |t| < 2 and non-critical for 
|r| > 2. This should affect the universal subleading power-law corrections to 
the asymptotic behavior. 

In [19], the following conjecture was made. For |r| < 2, the model corre- 
sponds to a theory with central charge c = r — 1 (corresponding to the analytic 
continuation of a model of r — 1 free bosons). This implies the following be- 
havior for the series 7p5^ counting coloured prime alternating tangles (with, 
say, four external legs): 



p^oo (^)_2 , . T -2- y/(2 - r)(26 - r) , , ^ ^ 
7p(r) ~ est gc[T) Pp^^ > 7(t) = 1"^' 

This was tested numerically in [19], but the results are not entirely conclusive 
(see also [10, 27]). 

In particular, as a corollary of the conjecture above, one would have the 
following asymptotic behavior for the number of prime alternating knots: 

/p ~ CSt5c(0) 6 

It is most likely that one can remove the "prime" property without changing 
the form of the asymptotic behavior (only the non-universal coefficient of the 
exponential growth would be modified); one can speculate that removing the 
"alternating" property will not change it either. 

Note the similarity between our problem of counting knots with that of 
counting meanders [C]. There too, the problem can be rewritten as a matrix 
model and the asymptotic behavior is dictated by 2D quantum gravity, leading 
to a non-rational critical exponent. A key difference is that in the case of 
meanders, corrections to the leading behaviour are expected to be power-law, 
making numerical checks reasonably easy. In contrast, if the conjecture above 
for knots is correct, the corrections are expected to be logarithmic (the theory 
being asymptotically free in the infra-red), which would make numerical checks 
extremely hard. 
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